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Abstract 

We consider a wireless relay network with one source, one relay and one destination, where 
communications between nodes are preformed via N orthogonal channels. This, for example, is the 
case when orthogonal frequency division multiplexing is employed for data communications. Since the 
power available at the source and relay is limited, we study optimal power allocation strategies at the 
source and relay in order to maximize the overall source-destination capacity under individual power 
constraints at the source and/or the relay. Depending on the availability of the channel state information 
at the source and rely, optimal power allocation strategies are performed at both the source and relay 
or only at the relay. Considering different setups for the problem, various optimization problems are 
formulated and solved. Some properties of the optimal solution are also proved. 



Index Terms 

Optimal power allocation strategy, amplify-and-forward relay networks, greedy algorithm, orthog- 
onal channels, sum capacity. 



I. Introduction 

Driven by the ever increasing demand for high-speed services within limited spectrum re- 
sources, cooperative relay transmission, as a promising spectrally efficient technique, has received 
significant interests in recent years [[Il-[|4l. Cooperative relaying promises numerous gains for 
wireless networks such as improved reliability [2] and increased network capacity [[31. The 
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benefits of relay transmission can also be reaped by orthogonal frequency division multiple 
(OFDM) access systems to support broadband services, e.g., the IEEE 802. 16j flU. 

Power limitation is a common problem in wireless networks. Therefore, allocating the limited 
power resources to the network nodes (i.e., the source and relay nodes) is a design consideration 
which has received much attention Il5ll- [|2T]| . For example, relay selection schemes, where only 
a subset of relay nodes is considered for cooperation, are studied in [I5]|-[|9l| as a simple power 
allocation technique. It has been shown in [ITOll that the optimal power allocation between the 
source and relay nodes can improve the network performance. Moreover, considering a multi- 
hop network in IfTTTl . the optimal allocation of power among the hops is shown to significantly 
improve the performance. Also, in IIT2ll - [|T5l . assuming that a common message is sent from the 
source node to multiple amplify-and-forward (AF) or decode-and-forward (DF) relay nodes 
over orthogonal channels, optimal power allocation strategies among relay nodes are studied. 

Recently, considering relay networks with multiple source-destination pairs (multicast), rele- 
vant power allocation strategies have been addressed lfT6ll - [l2TII . The set up of |fT6l is that one 
source sends independent messages to a number of destinations. A single OFDM channel is 
used between the source and a single relay node, where each private message is assigned to 
one frequency tone and then amplified and forwarded on different frequency tones to various 
users. Thus, power allocation at the relay is studied to maximize the minimum SNR among 
all destinations. In lfT7l - lfT9l , various power allocation schemes at both the source and relay 
nodes are developed for the general multi-source, -relay, and -destination set up. Transmissions 
are done over orthogonal channels, and the corresponding power allocation schemes are based 
on maximizing the minimum SNR among all users, minimizing the maximum transmit power 
over all sources, as well as maximizing the network throughput, the minimum rate among all 
users, and the weighted-sum of rates. Power allocation at the relay nodes only in the multi-relay 
multi-destination set up is considered also in [|20ll . and is based on maximizing the sum capacity 
of the network. In ll2T1l . power allocation at both the source and relay nodes in a set up similar 
to [[TSl . |fT9l is studied for maximizing the sum network capacity under power constraints for 
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orthogonal subchannels. 

In this paper, we consider a two-hop AF relay network, which consists of a single source 
node, a single destination, and a single relay nodes. We assume that the communications between 
nodes occur across orthogonal subchannels, where each subchannel is assigned independent 
information. Thus, the network capacity is the sum of the individual capacities of subchannels. 
This setup may correspond to the case, for example, when OFDM signalling is employed between 
the nodes. Our motivation for considering this setup is the widespread integration of OFDM into 
various wireless standards. Therefore, power allocation among the orthogonal channels (across 
frequency tones) is an important issue. 

Taking into account individual power constraints at the source and relay nodes, we study 
strategies for optimal allocation of the limited power among the orthogonal subchannels. 
Depending on the availability of channel state information (CSI) at the source and/or the relay, 
the power allocation is done at both the source node and the relay node or only at the relay 
node. The goal is to maximize the overall data rate of the network. 

We consider two cases for data forwarding at the relay. In one case, we assume the information 
received on the i-th source-relay subchannel is amplified and forwarded on the z-th relay- 
destination subchannel. In the other case, similar to (221, we allow the relay to switch the source 
message received on one source-relay subchannel to another relay-destination subchannel. With 
optimal power allocation, this strategy significantly improves the overall achievable rate. 

Power allocation for a similar OFDM-based relay network has been studied in ll23l . The main 
differences of our work are as follows, (i) With global CSI available at both the source and relay, 
we first prove that the optimal solution should satisfy a certain symmetry. This symmetry, in 
turn, allows us to directly find the optimal power allocation for both the source and relay. In 
contrast, in ll23l . an iterative optimization approach is studied. An iterative solution may not find 
the global optima in a reasonable number of iterations. Our approach, however, finds the optimal 
solution in one shot, (ii) We also consider the situation when only the relay has global CSI. In 
other words, the source only knows the source-relay channels. This is a more practical setup 
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as the source need not know the relay-destination channels. For this case, we propose a simple 
greedy algorithm to maximize the overall network achievable rate. It is also worth noting here 
that the optimization method of [23] needs global CSI and cannot be used for the latter case. 

Note that availability of global CSI (and therefore optimal power allocation) only at the 
relay has been considered in, for example, [|20l . However, the existing optimal power allocation 
schemes at the relay mostly focus on maximizing the sum capacity and it is not clear how the 
source (who lacks global CSI) should allocate rates across subchannels. Instead, our approach 
maximizes the achievable sum rate, while considering the outage. For this, we study how the 
source chooses its actual data rates across subchannels to avoid outage. Interestingly, at high 
SNR, the achievable rate is very close to that obtained in the case when global CSI is available 
to both source and relay. 

The rest of this paper is organized as follows. In Section |lll we outline the system model. Our 
optimization problems for maximizing the achievable sum rate are also defined. In Section [nil 
we propose optimal power allocation strategies at both the source and relay or only at the relay 
for different cases with respect to the availability of CSI. Some of the properties of the optimal 
solution are also studied. Section |IV] provides asymptotic analysis for optimal power allocation 
at high SNR. In Section |Vl numerical results are provided to demonstrate the benefits of the 
optimal power allocation. The paper is concluded in Section |VIl 

II. System model and problem definitions 
A. System and channel models 

Consider a two-hop relay network which consists of one source node {S), one destination node 
{D), and one relay node {R) with no line of sight between S and D. Thus, R is employed in order 
to assist the communications between S and D. Moreover, it is supposed that communication 
between nodes is done accross N orthogonal links. 

Two cases are considered for data forwarding at R. In the first case, the message received on 
the i-ih. S to R subchannel (SRi) is simply amplified and forwarded by the relay on the i-th R to 
D subchannel {RDi). This case is referred to as AF relaying. In the second case, as in [|22l| . we 
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assume that the relay sorts S to R subchannels and R to D subchannels based on their quality 
and assigns data on the i-th best S to R subchannel to the i-th best R to D subchannel. This 
amplify-sort-and-forward setup is referred to hereafter as ASF. For both AF and ASF cases, 
complete transmissions from S to D are composed of two phases: (i) S sends the messages 
while R listens; (ii) R forwards the message to D. 

Considering that each node is equipped with a single antenna, the received signal at R through 
SRi subchannel can be expressed as 

Vi = gi^/P^.mi + Hi (1) 

where rrii is the source message with a unit energy transmitted by S on the link SRi, Ps^ is the 
transmit power on the link SRi, gi is the associated flat Rayleigh fading channel coefficient for 
this link, which is a complex Gaussian random variable with zero mean and variance a^^, i.e., 
Qi ~ CAf{0,agJ, and Ui stands for a complex valued additive white Gaussian noise (AWGN) 
such that rii ~ CAf{0, 1). 

For AF relaying, at the end of the second phase, the received signal at D via subchannel RDi 
can be written as 

Xi = hi^/P^^^^== + n'. (2) 

where E{-} denotes the expected operation, Pr^ is the transmit power assigned to the link RDi, 
hi ~ CA/'(0, al^) is the channel coefficient for this link, and n[ ~ CN{Q, 1) is the AWGN. 

Moreover, for the ASF relay network, the received signal at D from R via the i-th. ordered 
link (denoted as RD(^i)) can be represented as 

X(i) = h(i\xl Pr... — , = + n',--. (3) 

where the subscript (z) stands for the z-th ordered link when sorting is based on \hi\ in the 
decreasing order such that |/i(i)| > ■ ■ ■ > \h{N)\- Also, Pr^^^ is the transmit power on the link 
RD(i^, and n'^-^ ~ CA/'(0, 1) is the AWGN. Note that in © can be found from © when aU 
SRi links are ordered to SR{i) and again sorted based on \gi\ in the decreasing order such that 

1^(1)1 > ■•■ > |^(iV)|- 
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This system can be viewed as subsystems operating on separate orthogonal suchchannels. 
For the i-th subsystem, we can define a source, a relay, and a destination denoted by Si, Ri, 
and Di respectively, where the communication from Si to Di is done over SRi-RDi. This point 
of view reduces our system to N parallel conventional single-relay systems each operating on a 
single channel, which helps our optimization approach later. Using existing results for the SNR 
of conventional relay systems ||24]| . the overall SNR of the i-th subsystem in the AF case is 





2 


Ps^gi\ 




\hi\ 


2 + 1 



A = „ I „ . o , (4) 



Similarly, in the ASF case, the corresponding SNR is 

B. Problem formulations 

Depending on the availability of the CSI at the relay and source nodes two different setups 
are considered: (i) when global CSI is available at both the source node and relay node, i.e., 
both S and R have knowledge of both gi and /ij, Vi. (ii) only local CSI is available at the source 
node, i.e., S knows only gi,\/i while R, knows gi and /ij, Vi. Thus, in both cases, the transmitter 
has the knowledge of the channel. It is worth mentioning that for water-filling protocols, it is 
typical in standards deploying OFDM to provide the CSI to the transmitter [4]. 

As mentioned earlier, we consider two different types of relaying: AF and ASF. Thus, a 
total of four different cases with respect to type of relaying and availability of global CSI can 
be considered. For these four cases, we study optimal power allocation strategies across 
orthogonal links at both S and R. For given finite energy budgets at each node, our goal is to 
maximize the achievable rate of the network. As can be seen from Q, the SNR of the link z, 
i.e., SRi-RDi, and, thus the capacity of this link depend on both Ps^ and Pji^,yi. Therefore, 
we seek the optimal power levels Pg. and Pr. , Vz that maximize the sum rate of the network. 
Finding the optimal solution for the ASF relay network can also be motivated in a similar way 
from ©. 

The optimal power allocation problems in these four cases can be stated as follows. 

June 29, 2010 DRAFT 



7 

1 ) Case I: AF relaying, global CSI at both S and R: The goal is to optimally allocate power 
in both S and R so that the sum capacity of the network is maximized under the individual 
power constraints at S and R. Thus, we have 

TV 

max > Ci (6) 

' ' 1=1 

N N 

s.t. J2Ps.<Ps, <Pij (7) 

1=1 i=l 

where Ci denotes the capacity of the link SRi — RDi, Ps and Pr are the individual power 
budgets at S and R, respectively. 

2) Case II: AF relying, global CSI only at R: In this case, we study power allocation strategy 
only at R. This is because, in this case, S knows only gi, Vi and, thus, the associated optimal 
power allocation at S is independent of hi,\fi. This means that the optimal power allocation at 
S is the traditional water-filling rule. Thus, we focus only on optimal power allocation at R. 

In the first phase, S decides the allocation of its limited power Ps to N orthogonal links using 
water-filling. For the obtained power allocation, then, S knows the link capacity for SRi, Vz. 
Based on the known link capacity, S decides a data rate 7^ for each link SRi, Vz. 

In the second phase, within its limited power Pji, the relay amplifies and forwards the source 
data received via links. Particularly, the data transmitted on SRi should be amplified (and 
forwarded on RDi) by the relay so that the overall capacity of this link SRi-RDi is equal to 
or greater than 7^. Otherwise, channel outage occurs. It is important to notice that the relay, due 
to its limited power, may not be able to forward all source's data. In such a case, no power 
should be assigned to the links whose data is not forwarded. In addition, on the links RDi that 
the relay decides to forward the source data, it is supposed to put just enough power in order 
to avoid the channel outage. The goal is, therefore, to forward as much source data as possible. 
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To achieve this goal, we can formulate the following optimization problem 



N 



max 

' i=i 

N 



(8) 

i=l 

s.t. 5^P,j, <Pr (9) 

i=l 

a e {7^,0} (10) 

where 7^ is the data rate decided by S for the link SRi, and Ci denotes the capacity of the link 
SRi-RDi controlled by Pr^ . 

The constraint (flOl) . as discussed earlier, is taken into consideration in order to avoid outage 
on links whose data is forwarded as well as to avoid power wastage on links whose data is not 
forwarded. Also, it is worth stressing that there is no point in allocating extra power on a link to 
increase Ci beyond 7/. The choice of 7j will be discussed in Section Hill where the optimization 
problem is solved. 

3) Case III: ASF relaying, global CSI at both S and R: Similar to Case I, we aim at finding 
the optimal solution for the power allocation at both S and R. However, the difference from 
Case I is that the data transmitted on the i-th ordered link SR^i) is amplified and forwarded on 
the i-th ordered link RD(^i). Therefore, we can pose the optimization problem as 



N 



max 



E^'w (11) 



N N 
i=l i=l 

where C(j) denotes the capacity of the link S'-R(j) — RD^iy 

4) Case IV: ASF relying, global CSI only at R: In this case, similar to Case II, the optimization 
problem can be formulated as 

N 



max 



i=l 

S.t. C'we{7w,o} (14) 



i=l 



June 29, 2010 



DRAFT 



where 7(i) is the data rate decided by S for the link S'i?(i), and C(i) denotes the capacity of the 
link S'-R(i)--Ri5(j) controlled by Puf^^y 

III. Optimal rate solution 
A. Case I: AF relaying, global CSI at both S and R 

In this case, it is assumed that both S and R know gi and hi, Thus, using dH), the individual 
capacity of link i is Ci = log (1 + pi) . Then, the sum capacity of the network is 

N N , I ,2|, |2 \ 

Based on ([151) . the optimal power allocation problem ©-([T]) can be reformulated as 

max,,. I^iog (i + _ :r:::T::ri . . i (i6) 



N N 



s.t. < 1, J]A < 1 (17) 

1=1 1=1 

where denotes the ratio of Ps, to Ps, i.e., = PsjPs, and similarly /3j = PrJ Pr. 

Because of constraints on a and (3, solving (fT6l)-(fT7l) directly can be difficult. Therefore, we 
first find a necessary condition on and /3j that significantly helps solving (fT6l)-(fT7l). For this, 
we notice that the individual power constraints in (flTI) lead to the following necessary condition 
on the sum of the given individual powers at S and R 

N N N 

Y,Ps.+Pr^<Ps + Pr or ^ + /3,r < 1 + r (18) 

i=l 1=1 i=l 

where r is the ratio of Pr to Ps, i.e., r = Pr/Ps- While this condition is not sufficient, any 
conclusion received from enforcing (fTSi) is necessary for the optimal solution of (fT6l) - (fT7l) . In 
the next theorem, we derive a relationship between optimal values of «j and 

Theorem 1: When and hi, \/i are known to both S and R, the optimal power allocation at 
S and R satisfies 

^_ a:Psh? + i 

a* P:PR\hi\^ + l 
where a* and /3* denote the optimal values of ctj and /3i, respectively. 
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Proof: Let Pj = [P5. PrJ denote a vector with elements of transmit powers used on 
subchannel i. Then, it can be shown from (fT5l) that for a given Pj, the sum capacity of the 
network is a concave function of Pj since d'^C/d'^Pi < for all values of Pj G 'R?,\/i. 
Therefore, (fT6l)-(fT7l) is a convex optimization problem with respect to Pj, Vi 

We now consider the objective (fT6l ) along with the constraint (fTSl ). The objective (fT6l ) can be 
optimized by using a Lagrangian multiplier method. The associated Lagrange function can be 
written as 

where A stands for the Lagrange multiplier. 

Taking derivatives of (|20|) with respect to Oj and /3j, Vi, and equating them to zero, we obtain 
that the optimal values of aj and Vi can be expressed as: 



«i = I 10 — + 7^W/3»^^ /3»^^ + T (21) 



1Ps\gi\^ 2V * V A 



2PR|/ii|2 2\/ r|/ij|2 V ' A 
Inserting (|22l) into (|2T1) . it is revealed that 



a*(a*P5|^?.r + 1) = r/3;(ArPi?J/i.r + 1) (23) 

and the claim of Theorem [T] follows straightforwardly. ■ 
Theorem [T] reveals an interesting symmetry for balancing power among individual links. Based 

on this result, if the optimal a* is known, then the optimal /3* can be immediately found. 
Noticing that (fT9l) is a necessary condition under which (fT6l) is optimized under the constraint 

( fTTl ). we can provide a modified optimization problem only based on ttj. So, the complexity of 

solving the optimization problem is greatly reduced. Then, using (fTS]) and (fT9l) . the sum capacity 

of the network can be reformulated with respect to only aj as 

(24) 



^ H a,Ps|^7.P + (l + A(a,))/2y' 
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where A{ai) = ^1 + 4:aiPs\hi\'^{aiPs\gi\'^ + 1). Then the problem (fT6l)-(fT7]) can be reformu- 
lated as 



where constraint (1271 ) is equivalent to X]i=i A — ^- Notice that in this optimization problem, 
the objective function and the constraints are all convex functions of a^. Therefore, traditional 
efficient numerical convex optimization techniques can be applied in order to solve it. 

B. Case II: AF relying, global CSI only at R 

In this case, as addressed earlier, we assume that in the first phase, S decides a power level 
P5- and a data rate 7^ for the link SRi. The goal of the relay in the second phase is to amplify 
and forward as much source data received over N source-relay links as possible, subject to 
avoiding outage on subchannels that R decides to forward through. This means that the data 
received over some links may not be forwarded by the relay, due to its limited available power. 
The optimization problem in this case is the problem (f8])-(fT0l). 

For given psRi, we can represent the individual achievable rate of the link SR—RDi as 




(25) 




(26) 



i=l 




(27) 




(28) 



The sum rate of the network can therefore be expressed as 



N N , 




) 



(29) 



where /3j is zero for links RDi that are not chosen to be amplified. 
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Based on ([8]) and (|29l) . subject to avoiding outage on likes with /3j ^ 0, the optimization 
problem can be reformulated as 



N 





\hi 


2 


PSR, + A^rI 


hi\ 


2 + 1 



max 

s.t. < 1 (31) 

1=1 

ae{7.,0},V^. (32) 

Due to the discrete constraint (|32l ). convex optimization techniques cannot be applied in order 
to solve (|30l)-(l32l). When seeking the optimal solution, it may be needed to consider all possible 
power allocations. The number of different power allocations, however, increases exponentially 
with A^. Therefore, we propose a greedy algorithm which finds near optimal solutions (and in 
many cases the optimal solution). The idea is to allocate the limited available power in the most 
efficient way. Particularly, notice that on the relay-destination link RDi, there is a minimum (3i 
guaranteeing successful source-destination communication. This value of (3i, denoted by /3i, is 
considered as the cost of communication on this link. This is because /3j is proportional to the 
amount of power spent on this link, if the link is chosen to be amplified and forwarded. Solving 



log 1 1+ :i . J =7. 

the value of A can be found as 

^ (27. _ 1) 27./^^ 

-27.)P«|/ii|2 ^^"^ 

where 5i denotes the ratio of ji to the capacity of the link SRi. While various subchannels can 

have different 5i, in this work we assume 5, = 5, Vz that is the worst-case. In Section |Vl we 

observe that even with 6i = 5, Vz, the achievable sum rate converges at high SNR to the case of 

Global CSI at both S and R. 

Consuming the cost (3i on link i, in return, the data rate 7j from the source to the destination 

is obtained. Then, we can define r/j = 7i//3i as the efficiency of allocating power to link i. Notice 

that the larger the value of rji, the better the link i. This is because for the same cost in terms of 
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the power spending, larger data rate 7^ is achievable on links that have larger values of 77^. So, 
all links are sorted in the decreasing order of 77^. We assign power to the links that have better 
r]i until we run out of power (before ^ • (3i becomes greater than one.) 

The greedy algorithm is designed such that the limited power is spent on the best links in terms 
of the achievable rate. However, this solution may not be optimal because by assigning power 
according to this greedy algorithm, we may end up with some positive leftover power. In such 
cases, other strategies that result in zero leftover power, may turn out to achieve a slightly higher 
rate. Nonetheless, the greedy algorithm is optimal in the sense of using the limited power in the 
most efficient way and in most cases provides the optimal solution for sum rate maximization 
problem. In particular, as grows large (i.e., the effects of the leftover power is negligible), it 
almost always gives the optimal solution. 

Now, we modify this greedy algorithm to avoid any leftover power. The idea is to allow 
S E [0, 1] to be a function of channel statistics. We notice that 6 remains constant as long 
as channel statistics are constant. To remove the leftover power, we consider the following 
max — max sum rate optimization problem 

N 



i=l 





\h. 




PsR, + PtPnl 


hi\ 


2+1 



(34) 



N 

s.t. J]A<1 (35) 

i=l 

aG{7^,0}, Vz. (36) 

To solve this optimization problem, for any given value of 5 G [0, 1], we use the aforementioned 
greedy algorithm to find the maximum sum rate. Then, we find such a 6* that provides the 
maximum sum rate. Please notice that in practice channel statistics do not change fast, so for 
given channel statistics, the optimal values of 6 can be found off-line and tabulated for later 
use. As soon as S* is chosen by the network, a single run of our greedy algorithm is needed to 
find the optimal power allocation per channel realization. 
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C. Case III: ASF relaying, global CSI at both S and R 

The only difference from Case I is channel ordering. In other words, data transmitted on 
SR(^i) (the i-th ordered link) is amplified and forwarded on RD(^iy The ordering is in the 
decreasing order of the channel gains \gi\'^ and respectively. Therefore, using dS]), the 

individual capacity of link (?) can be written as 

C(,)=log(l+P(,)). (37) 

where p(j) denotes the SNR received via link S'-R(i)--RD(j). The overall capacity of the network 
can be defined as 

Using (l38l) . the power allocation optimization problem (fTT1)-(fT2l) can be represented by the 
following problem 

max y log 1 H — (39) 

Af N 

s.t. J]«(.)<1, J]Ao<1- (40) 

i=l i=l 

Notice that the optimization problem (I39l)-(l40l) is mathematically equivalent to the problem 
(fT6l)-(fT7l) in Case I. The only difference from (fT6l)-(fT7l) is that the indices of optimal powers 
a^.^ and in (|39|)-(|40|) correspond to the ordered links. Thus, in order to find the optimal 
solution to (|39|)-(|401). we follow the same approach as in Case I by replacing gi, hi, ai and /3j 
with g^),h^i^,a{i) and respectively. 

D. Case IV: ASF relying, global CSI only at R 

Again, the optimization problem in this case is mathematically identical to the one in Case II. 
The only modification needed is to replace gi,hi, l3i and 7^ with gi^i), hi^i),(5i^i) and 7(j), respectively. 
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IV. Asymptotic analysis with global CSI at both S and R at high SNR 

In this section, we study the optimization problem of Case I at high SNR, i.e., when Ps and 
Pr tend to infinity. Again, we assume Pr = rPs- Thus, for a given r, we study the case that 
Ps tends to infinity. Through this analysis, we discuss the behavior of the optimal solution at 
high SNR regime. Some interesting observations are discussed. 

Using dH) and (|24|) . when Ps increases, the overall SNR on link i can be given by 



lim Pi = lim 

Pg— >oo Pg— >-o<: 



aiPs\9i 


m-l + 2a^Ps\ 


hM\)/2 


aiPs\9i\ 


2 + (1 + 2a,Ps 


\hi\ 





(41) 



= aiifiPs 

where v^j = + |/ij|/|(7j|). Then, it can be obtained asymptotically from (|4T1) that the 

sum capacity of the network is expressed only in terms of ctj by 

N N 

lim C= lim y^C, = y^\og{am) + N\og{Ps). (42) 

Pg— !>00 Pg— !>00 — ^ 
i=l 1=1 

Based on (|25l) and (|42)) . the optimization problem in Case I can be rewritten at high SNR as 

^ f P \ 

max V log {ai^i) + iV log —— (43) 
Oi.Vi ^ V 1 + r / 

1=1 ^ ' 

N 

s.t. ^ fti < 1 (44) 



i=l 



N 



l^-'M^'- '''' 

1=1 I *l 



It is worth noticing that the objective function (1431) is a monotonic function of a^, Vi. Using the 
Lagrange multiplier method, optimal solution to (|43l)-(|45l) can be obtained. Then the Lagrange 
function can be formulated as 

N / \ / N 

\ I I — , Ifl. l \ 

(46) 



where Ai and A2 are Lagrange multipliers associated with the constraints. Differentiating (1461) 
with respect to ai, Vz and equating the results to zero, we obtain the following closed-form 
expression for the optimal value of 

a* = max <( 0, + ^M] } ■ (47) 
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In (1471) . the Lagrange multipliers Ai and A2 are found once for a given set of values gi,\/i 
and hi,\/i. Therefore, it can be observed that when A2 > 0, for two subchannels with hi = hj, 
the optimal power allocation solutions put more source power on the weaker subchannel. In 
other words, if gi < gj then a* > a*, and vice versa. This is the inverse of a more traditional 
water-filling scheme, where only one power constraint exists and the optimal solution puts more 
power on the stronger subchannel. Similarly, rewriting the above equations in terms of /3*, we 
have 



which means that the behavior of the optimal (3* is the inverse of water-filling for A2 > 0. 

Furthermore, from Theorem [H at high SNR, it is obtained that the optimal solution to the 
problem (|43l)-(|45l) must satisfy 



Notice that Psct* and Pr/3* represent the allocated power to links SRi and RDi, respectively. It 
can be concluded then that when l^fjl > \hi\ on link i, the optimal assigned power on link SRi 
must be less than the optimal assigned power on link RDi. 



Consider a two-hop single relay network where nodes communicate through G {4, 20} 
orthogonal channels. The source-relay and relay-destination channels are assumed to be inde- 
pendent Rayliegh flat fading with variance a'^. The total sum power assigned to these 2N channels 
at both the source and relay is P = P5 + Pr. The noise is assumed to be additive zero-mean 
white Gaussian with unit variance. Thus, the total transmit SNR, which is used in our figures, 
is defined as SNR= Pa'^. Notice that transmit SNRs at S and R can individually be expressed 
as Pscr'^ and Prct^, correspondingly. 

A. AF relaying 

Assuming that the AF strategy is employed at the relay node. Fig. [2] shows the achieved sum 
capacity versus the total transmit SNR for Case I when the global CSI is available at both the 
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source and relay nodes and the proposed optimal power allocation strategy is used. Here r = 1 
(i.e., Ps = Pr = P/2). It can be seen from the figure that the sum capacity increases for 
larger N. Intuitively, this is because for limited Ps and Pr, larger N provides more chances of 
obtaining better source-destination subchannels on which the power allocation can result in a 
better sum capacity. 

Fig. |3] depicts the results for Case II, when the global CSI is available only at the relay and, 
thus, the power allocation can only be made at the relay. The sum rate is shown versus the total 
transmit SNR for r = 1 and = 20. Three scenarios for the greedy algorithm are considered. In 
the first one, 5 = 0.5 is fixed for all SNRs, while in the second and third scenarios, the optimized 
S, i.e., 6*, is used. The difference between the second and third scenarios is that water-filling at 
the source is assumed in one, while equal power allocation across all subchannels at the source is 
used in the other. The proposed greedy algorithm is used in order to allocate Pr to orthogonal 
subchannels. It can be seen from the figure that, as expected, higher data rates are achieved for 
optimal 6*. Interestingly, the equal power allocation at the source does not suffer from significant 
rate loss. Thus, in further examples for greedy-based power allocation, we consider only equal 
power allocation at the source. 

In Fig. m the results of the optimal power allocation in Cases I and II are compared to each 
other in terms of the sum rate plotted versus the total transmit SNR when r = 1, = 20, and 
optimal 5* is used. In agreement with our analysis, the global optimization in Case I always 
outperforms the other case. Particularly, at sum rate of 2 bits/s/Hz, Case I is shown to obtain 
4.5 dB power gain compared to Case II. Here, it is worth mentioning that this significant power 
gain is achieved at a modest value of A^. 

The results of the optimal power allocation in Cases I and 11 are also compared to each other 
in Fig. [5] in terms of the sum rates when r = 0.5 (i.e., Ps = 2Pr), N = 20, and 5 = 5*. It 
can be seen from the figure that again. Case I outperforms Case II for all values of the total 
transmit SNR. However, the performance gap between Cases I and II decreases as compared 
to the set up examined in Fig. HI This is because for r ^ 1 all source-relay subchannels are 
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likely to have much better conditions than the relay-destination subchannels, while the sum rate 
is mainly affected by the power allocation across the relay-destination subchannels. In addition, 
in both figures (Figs. |4] and Fig. [S]), one can see that Case II converges to Case I at high SNRs. 
This is because all subchannels become almost deterministic at high SNR, and therefore, the 
global optimization of the power allocation is less beneficial, while optimization of the power 
allocation only at the relay used in Case II tends to be almost sufficient. This is because Case II 
takes into consideration as many best relay-destination subchannels as possible within a limited 
power budget. 

Let us now consider r = 2 (i.e., 2Ps = Pr). The corresponding sum rates is shown versus 
the total transmit SNR for both Cases I and II in Fig. [6]for = 20 and 6 = 6*. It can be seen 
from this figure that Case I is better than Case II for all values of SNR. Interestingly, comparing 
all cases tested, i.e., r = 0.5 (Fig. [5]), r = 1 (Fig. H]), and r = 2 (Fig. (6]), the global optimal 
power allocation is most beneficial when r = 1. 

B. ASF relaying 

Fig. |7] depicts the sum rate versus the total transmit SNR for Case III when the global CSI is 
available at both the source and relay nodes and the optimal power allocation is performed for 
ASF relaying. For this figure, N E {4, 20} and r = 1. Case I is also depicted for comparison 
purposes. It can be observed from this figure that Case III outperforms Case I, as expected. 

Optimizing 6 and the power allocation only at the relay in Case IV, the sum rate is shown in 
Fig. [8] versus the total transmit SNR for = 20 and r = 1. The sum rate curves for Cases I, 
II, and III are also depicted for comparison. It is confirmed in Fig. [8] that Case III always 
outperforms all other cases including Case IV for the same reasons that have been explained 
while comparing Cases I and II. Interestingly, it can also be observed in this figure that Case IV 
is superior to the global optimization of Case I for moderate and high SNRs. 
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VI. Conclusions 

Two-hop AF and ASF relay networks, consisting of single source, destination, and relay nodes 
with limited individual power constraints at the source and relay are considered. Optimal power 
allocation strategies across N orthogonal subchannels between nodes are studied, depending on 
the availability of CSI at both the source and relay or only at the relay. 

When source has global CSI knowledge, through a symmetry property proved for the optimal 
power allocation, the optimization problem is solved. Via an asymptotic analysis, it is found 
that the global power optimization assigns more power on the weaker subchannels, which is 
the inverse of the traditional water-filing. When only the relay has global CSI knowledge, a 
greedy algorithm maximizing the achievable sum rate is proposed. For this, the actual data rate 
at the source is optimized and the minimum powers on subchannels are found at the relay in 
order to guarantee successful source-destination communication. It is also shown numerically 
that the optimal power allocation performed only at the relay can outperform the global power 
optimization scheme at moderate and large SNRs if simple subchannel sorting capabilities are 
added at the relay. 
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//orthogonal channels //orthogonal channels 



Fig. 1. Block diagram of a two-liop single relay wireless network with A'^ source-relay and A'^ relay-destination orthogonal 
channels. 
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-^Case 1, N=4, optimal power allocation at both S and R 
-B-Case 1, N=20, optimal power allocation at both S and R 
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Fig. 2. The sum capacity versus the total transmit SNR for Case I: the optimal power allocation at both the source and the 
relay for N G {4, 20} and r = 1. 
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Fig. 3. The sum rate versus the total transmit SNR using the proposed greedy algorithm for = 20 and r = 1. Three cases 
of S are examined: (i) optimal S* with the water-filling at the source (WF-based optimal 5*); (ii) optimal S* with the equal 
power allocation across A'^ orthogonal channels at the source (EP-based optimal S*); (iii) fixed 5 = 0.5 with the equal power 
allocation across A'^ orthogonal channels at the source (EP-based fixed S = 0.5). 




Fig. 4. Comparison of the sum rates for Cases I and II: A'^ = 20, r = 1, and S = 5* for equal power allocation across A'^ 
orthogonal channels at the source. 
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-B-Case 1, optimal power allocation at both S and R 
-©"Case II, optimal 8 , greedy-based power allocation 
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Fig. 5. Comparison of thie sum rates for Cases I and II: = 20, r = 0.5, and h = 5* for equal power allocation across A'^ 
orthogonal channels at the source. 
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-H-Case 1, optimal power allocation at both S and R 
-©"Case II, optimal 6 , greedy-based power allocatior 
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Fig. 6. Comparison of the sum rates for Cases I and II: A*" = 20, r = 2, and (5 = (5* for equal power allocation across A'^ 
orthogonal channels at the source. 
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-*-Case III, N=4, optimal power allocation at both S and R 
-HCase III, N=20, optimal power allocation at botli S and R 
-0-Case 1, N=4, optimal power allocation at both S and R 
-B-Case 1, N=20, optimal power allocation at both S and R 
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Fig. 7. The sum capacity versus the total transmit SNR: Case III of the optimal power allocation at both the source and the 
relay for A'^ £ {4, 20}, and r = 1. Case I is depicted for comparison. 



-B-Case I, optimal power allocation at both S and R 
-^Case II, optimal 8 , greedy-based power allocation 
-^Case III, optimal power allocation at both S and R 
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Fig. 8. The sum rate versus the total transmit SNR using the proposed greedy algorithm in Case IV for N — 20, r = 1, 
S = 5* for equal power allocation across A*' orthogonal channels at the source. Cases I, II, and III are depicted for comparison. 
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